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Abstract. We consider the free boundary problem for current- vortex sheets in ideal incompressible 
magneto-hydrodynamics. It is known that current-vortex sheets may be at most weakly (neutrally) 
stable due to the existence of surface waves solutions to the linearized equations. The existence of such 
waves may yield a loss of derivatives in the energy estimate of the solution with respect to the source 
terms. However, under a suitable stability condition satisfied at each point of the initial discontinuity 
and a flatness condition on the initial front, we prove an a priori estimate in Sobolev spaces for smooth 
solutions with no loss of derivatives. The result of this paper gives some hope for proving the local 
existence of smooth current-vortex sheets without resorting to a Nash-Moser iteration. Such result 
would be a rigorous confirmation of the stabilizing effect of the magnetic field on Kelvin-Helmholtz 
instabilities, which is well known in astrophysics. 

1. Introduction 

1.1. The Eulerian description. We consider the equations of incompressible magneto- hydrodynamics 
(MHD), i.e. the equations governing the motion of a perfectly conducting inviscid incompressible plasma. 
In the case of a homogeneous plasma (the density p = const > 0), the equations in a dimensionless form 
read: 

' d t u + V • (u <g) u - H <8 H ) + V<7 = . 
d t H - V x (u x H ) = , 
divu = 0, divH = 0, 

where u = (u\, 142,1*3) denotes the plasma velocity, H = (Hi, H2, H3) is the magnetic field (in Alfvcn 
velocity units), q = p + \H\ 2 /2 is the total pressure, p being the pressure. 
For smooth solutions, system (1) can be written in equivalent form as 

'd t u + (u ■ V)u - (H ■ V)H + Vq = 0, 

d t H + (u ■ V)H - (H ■ V)u = , (2) 
divu = 0, dWH = 0. 

We are interested in weak solutions of (1) that are smooth on either side of a smooth hypersurface 
T(t) = {x-3 = f(t, x')} in [0, T] x f2, where Q C M 3 , x' = (xi, X2) and that satisfy suitable jump conditions 
at each point of the front T(i). For simplicity we assume that the density is the same constant on either 
side of T(t). 
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Let us denote S^J ± (£) = {x 3 ^ f(t,x')}, where fl = Sl + (t) U f2 (£)ur(f); given any function g we denote 
g ± = g in Q ± (t) and [g] = gi, — gZ, the jump across T(t). 

We look for smooth solutions (it*, i? ± , q^) of (2) in f2 ± (i) such that T(t) is a tangential discontinuity, 
namely the plasma does not flow through the discontinuity front and the magnetic field is tangent to 
r(f), see e.g. [8], so that the boundary conditions take the form 

a = u ± ■ n , H ± ■ n = , [q] = on T(t) . 

Here n = n{t) denotes the outward unit normal on dil~ (t) and a denotes the velocity of propagation of 
the interface T{£). With our parametrization of T(t), an equivalent formulation of these jump conditions 
is 

d t f = u ± -N, H ± -N = ) [q] = on T(t) , (3) 

with N := (—dif,—d2f,l). Notice that the function / describing the discontinuity front is part of the 
unknown of the problem, i.e. this is a free boundary problem. 
System (2), (3) is supplemented with initial conditions 

u±(0,x) = u±(x), H ± (0,x) = H±(x), zeO^O), 

f(o,x') = f (x>), z'er(o), w 

where divu^ = divH^ = in 0,^(0). The aim of this article is to show a priori estimates for smooth 
solutions to (2), (3), (4). This must be seen as a preliminary step before proving the existence and 
uniqueness of solutions to (2), (3), (4). 

In the last years there has been a renewed interest for the analysis of free interface problems in fluid 
dynamics, especially for the Euler equations in vacuum and the water waves problem, see [6, 7] and the 
references thereinto. This fact has produced different methodologies for obtaining a priori estimates and 
the proof of existence of solutions. If the interface moves with the velocity of fluid particles, a natural 
approach consists in the introduction of Lagrangian coordinates, that reduces the original problem to a 
new one on a fixed domain. This approach has been recently employed with success in a series of papers 
by Coutand and Shkoller on the incompressible and compressible Euler equations in vacuum, see [6, 7]. 
However, this method seems hardly applicable to problem (2), (3), (4). 

In the present paper we follow a different approach. To reduce our free boundary problem to the fixed 
domain, we consider a change of variables inspired from Lannes [9] . The control of the function describing 
the free interface follows from a stability condition introduced by Trakhinin in [13]. The a priori estimate 
in Sobolev norm of the solution is then obtained by showing the boundedness of a higher-order energy 
functional. 

1.2. The reference domain ft. To avoid using local coordinate charts necessary for arbitrary geome- 
tries, and for simplicity, we will assume that the space domain f2 occupied by the fluid is given by 

n := {(x 1: x 2 ,x 3 ) Gl 3 |i'= {x 1 ,x 2 ) e T 2 , x 3 e (-1, 1)} , 

where T 2 denotes the 2-torus, which can be thougt of as the unit square with periodic boundary conditions. 
This permits the use of one global Cartesian coordinates system. We also set 

fi ± := fi n {x 3 ^ 0} , r := n n {x 3 = 0} . 

On the top and bottom boundaries 

r± :={(x',±l),*'GT 2 } 
of the domain fi, we prescribe the usual boundary conditions 

u 3 = H 3 = on[0,T]xr±. (5) 
The moving discontinuity front is given by 

r(t) := {(x', x 3 ) £T 2 xI,i 3 = fit, x 1 )} , 



where it is assumed that —1 < f(t, •) < 1. 
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1.3. An equivalent formulation in the fixed domain ft. To reduce the free boundary problem (2), 
(3), (4), (5) to the fixed domain 0, we introduce a suitable change of variables that is inspired from [9]. 
This choice is motivated below. In all what follows, H s (u>) denotes the Sobolev space of order s on a 
domain uj. We recall that on the torus T 2 , H S (T 2 ) can be defined by means of the Fourier coefficients 
and coincides with the set of distributions u such that 



J2 (i + mtic„hi 2 <+^. 



c n (u) denoting the n-th Fourier coefficient of u. The following Lemma shows how to lift functions from 

r to n. 

Lemma 1 ([9]). Let in > 1 be an integer. Then there exists a continuous linear map f € -ff™~°' 5 (F) i— > 
ijj € H m (il) such thati/)(x',0) = f(x') on T, ip(x',±l) = on T ± , and moreover d 3 tp(x' , 0) = 0i/m>2. 

For the sake of completeness, we recall the proof of Lemma 1 in Section 7 at the end of this article. The 
following Lemma gives the time-dependent version of Lemma 1. 

Lemma 2. Let m > 1 be an integer and let T > 0. Then there exists a continuous linear map 

f e n7 = -VJ([0,T];i/ m ^-°- 5 (T 2 )) ^ </< e n™" V J ([ . r ]^ m ^(^)) such that = f(t,x'), 

ij}(t,x' ,±1) = 0, and moreover d 3 ip(t,x' ,0) = if m > 2. Furthermore, there exists a constant C > 
that is independent of T and only depends on m, such that 

V/en^V J '([o,T];F m -^- - 5 (T 2 )), Vj = 0,...,ro-1, Vte[o,T], 

\\d}iit, -)\\ H m- jm < c \\dif(t, .)||k— i-o..(r») ■ 

The proof of Lemma 2 is also postponed to Section 7. The diffcomorphism that reduces the free boundary 
problem (2), (3), (4), (5) to the fixed domain f2 is given in the following Lemma. 

Lemma 3. Let m > 3 be an integer. Then there exists a numerical constant Eq > such that for all 
T > 0, for all f e nf^ 1 ^ ([0,T]; H m ^-°- 5 (T 2 )) satisfying \\f\\ n[0 ,T]-m-^)) < e , the function 

*(t,x):=(x',x 3 +il>(t,x)), (t,x) e [0,T] x fi, (6) 

with ip as in Lemma 2, defines an H m -diffeomorphism of £1 for all t € [O,? 1 ]. Moreover, there holds 
di* G V([0,T];H m -i{n)) for j = 0,...,m-l 7 *(*,x',0) = (x', f{t, x')), *{t,a> ,±l) = (x',±l), 
d 3 <Z>(t,x',0) = (0,0,1), and 

V*e[0,T], ||V(V)lki.~(n) < \- 

Proof of Lemma 3. The proof follows directly from Lemma 2 and the Sobolev imbedding Theorem, be- 
cause 

d 3 V 3 (t,x) = l+d 3 1p(t,x) > 1- ||V(V)lk([0,T];Vl'i.~(fi)) > 1 - C||/|k([0,T];H^(T2)) > 1/2, 

provided that / is taken sufficiently small in ^([0, T]; H 2 - 5 (T 2 )). In the latter inequality, C denotes a 
numerical constant. The other properties of ^ follow directly from Lemma 2. □ 



We set 



and we compute 



A := [D3 r ] _1 (inverse of the Jacobian matrix) , 

J := dct [D^] (determinant of the Jacobian matrix) 

a := J A (transpose of the cofactor matrix) , 





-9nl> 
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(9) 



We already observe that under the smallness condition of Lemma 3, all coordinates of A are bounded by 
2 and J £ [1/2; 3/2]. Now we may reduce the free boundary problem (2), (3), (4), (5) to a problem in 
the fixed domain ft by the change of variables (6). Let us set 

v ± (t,x) :=u ± (t,^{t,x)), B ± (t,x) := Q ± (t,x) := q ± (t, xj) . 

Then system (2), (3), (4), (5) can be reformulated on the fixed reference domain f2 as 

'dtv* + (w ± • V)w ± - (B± • V)B ± + A T VQ± = , 
d t B ± + (v* ■ V)B ± - (B ± ■ V)« ± = , 

(A T V)-V ± = 0, (A T V)-B ± =0, m[0,T}xn ± , 

d t f = v ± -N, B ± -N = 0, [Q} = 0, on[0,T]xL, (8) 

vf=B^=0, on[0,T]xL±, 

J\t=o = fo, onT. 
In (8), we have set 

JV := (-5^,-92^,1), 

v := Av — (0, 0, W/J) = (»!, v 2 , (v-N- d t ip)/J) , B := A B = (B u B 2 , B ■ N/J) . 

Vectors are written indifferently in rows or columns in order to simplify the redaction. Notice that 

J = l, N = (-dif, -d 2 f, 1) onT, £3=^3 = onLandr±. (10) 

We warn the reader that in (8), the notation A T is used to denote the transpose of A and has nothing 
to do with the time interval [0, T] on which the smooth solution is sought. We hope that this does not 
create any confusion. 

1.4. The main result. 

1.4.1. The linearized stability conditions. The necessary and sufficient linear stability conditions for pla- 
nar (constant coefficients) current-vortex sheets was found a long time ago by Syrovatskii [12] and Axford 
[2]. Let us consider constant vectors w, ± ,i? ± satisfying (3) with the planar front f(t,x') = at + £' • x 1 
and constant pressures q ± = 0. (Here we consider for this paragraph that x' belongs to K 2 instead of T 2 
and X3 G R. This is however of no consequence on what follows.) The linear stability conditions for such 
piecewise constant solutions to (1) read 

N| 2 <2(V+| 2 + |!n 2 ), (lla) 

\H+ x [u]\ 2 + \H- x [u]\ 2 < 2 \H+ x H-\ 2 . (lib) 

Under the additional assumption H + x H~ ^ 0, then (lla) follows from (lib) and the strict inequality in 
(lla) follows from the strict inequality in (lib). The case of equality in (lib) corresponds to the transition 
to violent instability, i.e. ill-poscdncss of the linearized problem. In the region of parameters defined by 
(11). the associated linearized equations admit surface waves of the form exp(i rt + in ■ x' — \n\ [0^3 [) for 
n £ M 2 \ {0} and some suitable r £ R, see [12, 2] or [4, page 510]. We also refer to [1] for the derivation 
of weakly nonlinear surface waves. 

The interior of the set of parameters described by (11) is defined by the condition 

\H+ x [u]\ 2 + \H~ x [u]\ 2 < 2 \H+ x H-\ 2 . (12) 

In particular, H + x H~ 7^ and (lla) becomes redundant. The condition (12) is always satisfied for 
current sheets, i.e. if [u] = and H + x H ~ 7^ 0. If [u] 7^ 0, condition (12) can be rewritten as 

V2\H+\\H-\\ S m(<P + -<p-)\ 

I N I < / , 

^J\H+\ 2 sin 2 (p+ + \H-\ 2 sin 2 99- 
where ip^ denotes the oriented angle between [u] and H ± . 
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Under the "spectral stability condition" (12), Morando, Trakhinin and Trebeschi [10] have shown an 
a priori estimate with a loss of three derivatives for solutions to the linearized equations with constant 
coefficients. In this paper we shall consider the following more restrictive situation: 

max(\H + x [u]\, \H~ x [tt] Q < \H + x H~\. (13) 

Under the latter more restrictive stability condition, which represents "half" of the stability domain 
defined by (12), Trakhinin [11] has shown an a priori estimate in the anisotropic space Hi, without loss of 
derivatives from the data, for solutions of the linearized incompressible equations with variable coefficients. 
Similar stability conditions have also been considered by Trakhinin for the analysis of linearized and 
nonlinear stability of compressible current-vortex sheets, see [13, 15, 5]. The choice of the space Hi in 
[11] was motivated by the fact that the free boundary T(t) is characteristic. However, we shall prove here 
that no loss of derivatives in the normal direction to the boundary occurs and we shall obtain estimates 
in standard Sobolcv spaces. Though there is no loss of derivatives from the source terms of the equations 
to the solution in the main a priori estimate of [14], the regularity assumptions on the coefficients are 
rather strong (stronger than what we shall assume here), and it is not so clear that the estimate in Hi 
is sufficient to prove an estimate in some H™, m large enough, with coefficients in the same space H™. 
There are even strong reasons to believe that with the formulation of [14], a loss of regularity will occur 
with respect to the coefficients of the linearized equations. 

Our goal here is to prove a closed estimate where coefficients are estimated in the same space as the 
data. As a matter of fact, we have found it more convenient to work directly on solutions to the nonlinear 
equations. Since we are considering classical solutions in three space dimensions, our a priori estimate 
will be proved in H 3 (£l), a space that is imbedded in W 1,oa by the Sobolev imbedding Theorem. 

1.4.2. The main result. For a pair of functions u = (u + ,u~) € £P(S! + ) x iJ s (Q _ ), with real s > 1, we 
will shortly write 

||u + |U+ := ||« + ||^+) , IMk- := |M|^ (n -) , H^lki := \\u + \\ s , + + |Mk- • 

We also let | • k± denote the LP norm on Vt^, and | ■ \ p denote the LP norm on il for p > 1 and p ^ 2; 
the L 2 norm on f2 ± is denoted by || • Ik. Our main result reads as follows. 

Theorem 4. Let S Q e]0, 1/2], let R>0, and let v^,B^ e i? 4 (fi ± ), f e H 4 - 5 (T 2 ) satisfy 

Vz'eT 2 , \B+xBo(x',0)\>5 , 
\fx' e T 2 , max(|B+ x [v ] (x',0)\, \B„ x [v ] {x',0)\) < (1 - So) \B+ x B " {x',0)\ , (14) 

ll^lki + ll^lki + II/o||h3.5 (t2) < r. 

Then there exist e± > 0, To > and C\ > that depend only on 5q and R such that if \\fo\\H 2 ^n 2 ) ^ £ i> 
then for all solution (v ± ,B ± ,Q ± ) g <g([0, T]; H 4 ^)), f e %f([0, T}; H 4 - 5 {T 2 )) to (8) satisfying (without 
loss of generality) 

/ Q~(t,x)dx+ Q + (t,x)dx = 0, 
Jo.- Jn+ 
for all t € [0,T], the following estimates hold: 

ll^WIki + \\B±{t)h,± + ||Q ± (t)!l3,±+ \\f(t)\\ H 3., m < Ci , 

\\f(f)\\H"<J>) <2ei, 1 j 

for all t € [0, min{T, To}] . 

Directly from (8) and (15) it readily follows a uniform estimate for ||9tf ± (i)||2,±: ||5ti? ± (t)||2,± and 
\\dtf(t)\\ H 2.5 {T2) . 

The first two conditions (14) are nothing but a uniform version of (13) on the initial front. Then 
our main result gives a uniform control of solutions to (8) provided that a flatness condition is satisfied 
by the initial front. The main result also shows that the front remains sufficiently flat on a small time 
interval. The main interest of Theorem 4 is to show that energy estimates without loss of derivatives 
can be proved for (8) in the framework of standard Sobolev spaces. We hope that in a near future, our 
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approach will yield an existence and uniqueness result for (8) without using a Nash-Moscr iteration. As 
far as we know, no existence result has been proved yet for (8), with or without a Nash-Moscr iteration. 

1.4.3. Strategy of the proof. We consider the following energy functional 

*{t) := \\v ± (t),B ± (t)\\l ± + ||Q ± (t)||i >± + ||/(t)||^. 5(T2) + \\d t f(t)\\ 2 H2 . HT2) . (16) 

Even though this function is not conserved, it is possible to show that sup tg [ To ] ${t) remains uniformly 
bounded for sufficiently smooth solutions to (8), whenever To > is taken sufficiently small (To being 
independent of the solution that we are considering). The strategy for proving Theorem 4 is the follow- 
ing: we first estimate the velocity and magnetic field by showing energy estimates on their tangential 
derivatives (meaning the d\ and 82 derivatives), on their divergence and on their curl. Computing the 
curl equation is the crucial point if one wants to use standard Sobolev spaces (this is one difference with 
[14]). The front / will be estimated directly from the boundary conditions in (8). Eventually, the pressure 
will be estimated by showing that satisfy an elliptic system with source terms depending only on 
, B^ , / which have been estimated previously. Then we shall combine all these estimates to show that 
they yield a uniform control of solutions on a time interval that only depends on the size of the initial 
data. 

Not so surprisingly, Theorem 4 requires an additional degree of regularity on the solution compared to 
the space in which we prove the estimate. This technical point is assumed only to justify all computations 
below (integration by parts and so on). This is exactly the same as when one proves a priori estimates 
for solutions to first order hyperbolic problems and in many aspects our analysis is closely linked to 
techniques used in hyperbolic boundary problems with characteristic boundaries. In particular, if we 
believe that coefficients of the differential operators in (8) should have the same regularity as the solution 
to (8), then A should belong to H 3 if v ± ,B ± belong to H 3 . This forces the lifting tp of the front / to 
belong to H 4 and this is where it is crucial to gain half-derivative from / to if). This is the reason why 
we have adopted the same lifting procedure as in [9]. 

2. Estimate of tangential derivatives 

2.1. Uniform control of low order derivatives. From now on we consider a time T' > such that 
we have for our given solution the uniform estimates: 

Vte[0,T'], ||/(V)||if=5(T=) <eo, (17a) 

Wv^t) - ttf, B ± (t) - B±|| a>± < £o , (17b) 

where in (17), the numerical constant Eq is given by Lemma 3. Let us already observe that with our 
choice of £0, (17a) implies 

V(t,x)e[0,T']xn, |VV(*,ar)I < |. 

Moreover, the Sobolev imbedding Theorem implies that the H norm dominates the L°° norm on 
so we can further restrict £0, depending only on So, such that the following inequalities are implied by 
(17b): 

V (t, x') G [0, T'] x T 2 , \B+ x B~ (t,x',0)\ > -± , (18a) 

, , % r , max(|B+ x [v] (t,x',0)\, \B~ x [v] (t,x',0)\) S , , 

VforOeMxT*. — U ^4^)^ j ' j ^f- d8b) 

Of course, the time T' chosen above a priori depends on the particular solution that we are considering, 
and one of our goals is to show below that T' can be chosen to depend only on So and on the norm R of 
the initial data. 

We will denote generic numerical constants (for instance constants that appear in Sobolev imbeddings) 
by the same letter C or by Mq. Such constants are allowed to depend only on Sq and R. We also let 
F denote a generic nonncgative nondecreasing function which does not depend on the solution. In 
particular, we feel free to use F + F = F,FxF = F and so on. We shall sometimes write u(t) instead 
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of it(i, •), for some given function u depending on t and x. For shortness we shall write || ■u =1= , .B^ |[ 3,± 
for ||v ± ||3.± + ||B ± ||3 i ±, and similarly for ||9tf ± , St-B^ltai an d other quantities. Let us now turn to the 
derivation of L 2 estimates for tangential derivatives of the velocity and magnetic field. 

2.2. Estimates of tangential derivatives. Let us denote by 8 — (81,82) the horizontal (tangential) 
derivatives. Inspired from [13, 14] we define on [0, T] the energy functional 

^ ? eI(_v 

where X^ = X(v , B ± ) is a Iff 1 function that will be chosen appropriately later on. In particular, the 
choice of A ± will be made so that we have 

l^ + IU^([o,T']xr2+) < 1, P _ ||l^([o,t']xo-) < 1, (20) 

which will imply that the matrix in the integrals defining J4?(t) is positive definite (hence we shall recover 
a control of the tangential derivatives of the solution). 
We compute the time derivative 



± |a|<3 



e £ L (-a± • (tb±) dx 



± |q|<3 



E E J ± 8t\ ± 8 a v ± ■TB ± Ax 



± |q|<3 



± |a|<3 







(-;. 


-;•) ( 



{(S* • V)^ - (B ± • V)B ± + A T VQ ± }\ /g Q w ± 



EEL -A± 1 W,*± ™ D± ,6± ™„.±1 ' S°R± 



n± V- A± 1 / ^ 9 a {(^•V)B ± -(B ± -V) l . ± } j \S a B 



p=i 

where each term in the decomposition will be defined below, and we leave as a very simple exercise 
to the reader to check that the sum of all these terms coincides with the time derivative J^'(t). We now 
define and estimate all the terms in the decomposition of JK"(t). We first consider 



J^(t):=-Vy / 8 t X ± 8 a v ± ■8 a B ± dx. 



± \ct\<3' 

which is trivially estimated by 

ViG[0,T'], \Jf 1 (t)\<C^(t)Y / \\9tX ± \\ L ^(ni)- (22) 

± 

Next we consider some of the terms with the highest number of derivatives. Let us define 

This term is estimated by integrating by parts and recalling the boundary condition (10). We obtain 



^(*)=EE J i (divw ± +div(A ± S ± )) (|9 Q w ± | 2 + |9 Q B ± 

± |q|<3 ^ 

- (divi^ +div(A ± 5 ± )) 9 Q v ± ■ cf dx , 



s 



J.-F. COULOMBEL, A. MORANDO, P. SECCHI AND P. TREBESCHI 



from which we already get 

< CS{t) ||divt5 ± J divB ± || i « (n ±) + ||div (A ± 5 ± ), div (A ± B ± )|| L oo (n±) . 
± 

Using the expression of w ± , B ± , we get (recall that the estimate (17a) implies in particular l + d^ip > 1/2) 
Vie[0,T'], ||div6 ± ,divB ± || i » (n ± ) <F{g{t)), 

||div (A ± u ± ), div (A ± S ± )|| i « (n ±) < F(^(t)) ^±11^1.-^) ■ 

We thus end up with 

VtG[0,T'], |^ 2 (t)|<^(^(t))(l + ^||A ± || H , 1 , 00(0±) ). (23) 
Let us now consider the term 

= - E E / ATy ( r • { r w± - A± r B± } da; • 

This is the term which requires the most careful analysis. We first observe that the term in the sum 
which corresponds to a = (no tangential derivative) is estimated in an elementary way by Cauchy- 
Schwarz inequality, and admits an upper bound that is the same as in (23). We thus feel free to slightly 
modify the definition of J#§ and from now on we only consider the sum over the multi-indices a satisfying 
1 < | ct | < 3. A first integration by parts gives (here we use Einstein's convention over repeated indices) 

■*!(*) = E / A ^Q + {^V - ^ + d a B+) dx' 

1<M<3 

- E / A 3i TQ+ {Tvf-\+TB+}dx> 

i<M<3 r + 

- E [ A 3i TQ- [Tvr -\-TBr} dx' (24) 
i<M<3 r 

+ E / A 3l d a Q- {Tvr -\-d a Br} dx' 
i<H<3 r - 

+E E / ^Q^^A^rvt-x^Bt^dx. 

± 1<|q|<3 S1± 

Let us notice first that 

A 3l {Tvf - \±d a Bf}\ X3 = ±1 = j {Tvf - A ± r J B 3 ± }| X3=±1 =0, 

because of (7) and tp = v 3 = Bf = on [0, T] x T±. Therefore the second and fourth boundary integrals 
on r± in (24) vanish identically. As for the two boundary integrals on T, from (7), (10) and the boundary 
condition [Q] = on V we have 

A 3 . = N , [d a Q] = on T . 
Therefore we may rewrite (24) as J(f 3 (t) = J4f 3 i(t) + ffljait) with 

je 31 (t):= ]T fd a Q[(d a v-\d a B)-N]dx', (25) 

l<|a|<3 r 

m,i{t) := E E / d i (&* v ? - X± ^ B f)} dx > (26) 

± 1<|q|<3 n± 
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where [■] in (25) still denotes the jump across T, and Q denotes the common trace of on F. 

Let us first consider the term (t), which is where the choice of A ± is made. The boundary conditions 
[v ■ N] = B ± ■ N = on T yield W([v ■ N]) = d°'(B ± ■ N) = on T. Therefore we may write 



^3i(i) = - E d a Q[{d a -N]-v-\[d a ;N]-B 

1<|q|<3 r 



dx . 



where [d ; N] denotes the commutator between d and the multiplication by N . This commutator can 
be written as a sum of the form 

\$*;N]=TN+ *^NT~ P , 

l<|/3|<|a|-l 

where * denotes some harmless numerical coefficient. Let us assume for the time being that we can 
construct A ± on [0, T'} x F that satisfy 



A+ B+ - X- Bi = [vi] , 
A+ B+ - X- S 2 - = [v 2 ] , 



(27) 



so that [v' — XB'] = 0, where we have set v' := (vi,V2) and so on. Then the decomposition of the 
commutator reduces J^i(t) to 



•*5i(*)= EE*/ 



1<M<3 1<|/3|<|q|-1 



d a Qd"Vf-(d a ^v'-Xd a ^B') dx'. 



where we have set V := (<9i, d 2 ) (here the indices ± do not play any role so we feel free to omit them). 
We now recall the following classical product estimate. 

Lemma 5. The product mapping H°- 5 (T 2 ) x i/ 1 - 5 (T 2 ) — > H°- 5 (T 2 ), (/,<?) ' — ► f g is continuous. 

We can now estimate each term in the above decomposition of J^i(t). In the case \a\ — \/3\ = 1, we 
get (use Lemma 5 for the product estimate and the fact that ff 15 (T 2 ) is an algebra) 



)ds' 


< c 


Tq 













aPv'f-@*- f 'v'-\ar- p B 



< C \\VQ\\ H1 . 



(T) 



_ff°- 5 (T 2 ) 



a v — Xo B 



i?°-5(r) 



ff 15 (r) 



<F(,r(t))(i + ^||A ± || ffl 



■>(T) 



In the case \a\ — |/3| > 2, which only happens for \a\ = 3 and |/3| = 1, we have 



trQ8rVf-[d v'-Xd B 1 dx 



)dx' 


< c 


Tq 










ff-°- 5 (r) 



H 05 (T) 



< C \\VQ\\ Hl _ Hr) 



d P Vf 



O V — AO B 



j?o-5(r) 



<F(^))(l + ^||A ± || ff , 5 



B(T) 



Summing all the estimates, we have obtained 

vte[o,T'}, |^3i(i)|<^(^W)(i + Ell A± ll« 1 - 5 (n 



(28) 



provided that we can construct X^ that satisfy (27). Let us therefore turn to the construction of these 
functions. 

We first observe that the boundary conditions (10) give 

Bf = Bf & 1 f + B$ d 2 f , [v 3 ] = [v 1 ]d 1 f+ [v 2 ] d 2 f, on T , 
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so (27) is equivalent to the relation 

[v] = A+ B+ - \- B~ on T . 

Using the lower bound (18a) on the time interval [0, T'], we know that (27) is a Cramer system (otherwise, 
B + and B~ would be colinear). Hence A are uniquely determined on [0,T'] x T and have the same 
regularity as v ± , B ± on the boundary T. Moreover, the latter relations give 

' ^ ' ' ^ = \B+ x B-r ' ' ' 2~' 

where we have used (18b). As in [13, 14], we extend A* to the domains f2 ± as functions that do not 
depend on the normal variable X3. Using time or tangential differentiation on the system (27), we can 
easily obtain the estimates 

Vie [0,2"], ||A ± || ff i. 5(r) + ||A ± || w i.« (n ±) + 11^11^^) <F{S{t)), (29) 

60 
2 ' 

The latter estimates on A ± simplify (22), (23) and (28), and give 

Vie [0,2"], \J%(t)\ + \je 2 (t)\+\J%i(t)\<F(£(t)). (30) 

We emphasize that in the estimate (30), the nondecreasing function F depends on Sq because the estimates 
on A ± depend on 5q, but F does not depend on the particular solution that we are considering. 

Let us now consider the term ^32 (i) in (26). We decompose ^32 (i) as ^32 (i) — ^321 (i) + ^322 (i), 
with 



Il^ ± lk° c (f2±) < 1 



*ki(t):=Y, E / d a Q ± {d j A ji )(d a v?-\ ± d a B?)dx, 

± l<|o|<3 n± 

^32 2 (*):=E E / d a Q ± A Jl d J {d a vf-\ ± d a Bf)Ax. 



± 1<H<3' 

The first term J^s2l{t) is easily estimated by applying Cauchy-Schwarz inequality and by using the L°° 
estimate of X^, see (29): 

Vie [0,7"], |-as 2 i(t)l<*W))- (31) 

As for ^322 (t), since we have the divergence constraint Aji djvf = A 3 i djBf = 0, we may write 

^32 2 (t) = - E E / f A ^ dAv t + A H ( d J X± )d a Bf - X ± [T; A 3l d 3 \Bf\ dx , 

± i<H<3' /n± 

where [•; •] still denotes the commutator. The latter terms are now estimated in a somehow brutal way 
by applying Cauchy-Schwarz inequality We recall that the H A norm of ip is controlled by the H 3 5 norm 
of / thanks to Lemma 1, and that commutators in L 2 are controlled by standard estimates which may 
be found for instance in [3, page 295]. Eventually we obtain 

Vie [0,7"], |^322W| < F{S{t)) . 

Combining with (31), and (30), we end up with 

Vie [0,7"], \jn(t)\ + \J^(t)\ + \J^(t)\<F(^(t)). (32) 

Going on with the estimate of the terms in the decomposition (21) of Jff"(t), we finally consider 

and 

-EE / [d a ;A T v]Q ± -{Tv ± -\ ± d a B ± } dx. 
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Indeed the reader can check that the relation (21) holds with the above definitions of . . . ,^5. Ap- 
plying again the classical commutator estimates and using once again the L°° estimates of A ± , we have 

Vte[0,T"], \JtH(t)\ + \J%(t)\<F(*(t)). (33) 

Combining (32) and (33), we have therefore derived the inequality 

VtG[0,T / ] > \Jf{t)\ <F(£(t)), 

for a given nonnegativc nondecreasing function F that is independent of the solution. Integrating from 
to t £ [0, T'\ and using the L°° bounds on A ± , we have already proved our main a priori estimate for 
tangential derivatives: 

Vte[0,T'], V \\Tv ± (t),TB ± (t)\\l < M + tF(max £{s)), (34) 

^ — ' M± 0<s<t 

\a\<3 

where Mq is a numerical constant that only depends on Sq and R (here we have used (29) to derive a 
lower bound for the positive definite matrix appearing in the definition of the energy functional Jff). 

3. Divergence and curl estimates for v and B 

3.1. Estimates for the divergence. In this section we derive suitable estimates for the divergence of 
w ± ,_B ± in Expanding the divergence constraint for v ± 1 we find that for each t £ [0, T"], there holds 

divf — d 3 v^ + d 2 vf — d 3 vz + - d 3 vf = 111 , 

from which the identity 

± Vt(j-d 3 v ± . ± 

divu = m S2 



readily follows. Since H (fi ) is an algebra, we get 



Vte[0,T'], ||div« ± (t)|| 3i ± < C 



||%« ± (*)||2,±<C||/(t)|| Ha .. (Ta) ||i; ± (t)||3 1 ± 



2 

± 



The analogue estimate for the divergence of B is obtained by following the same lines, and we have 
thus proved the a priori estimate 

V* G [0,T'] , ||div» ± (t),divfl ± (t)|| 3l ± < Co ||/(i)|| ff 2.5 (T2) ||i; ± (i),B ± (t)||3 1 ± ■ (35) 

3.2. Estimates for the curl. In order to estimate the curl of v , .B we proceed as follows. Let us 
introduce the curl of the Eulcrian velocity and magnetic fields u, H 

£ := curl it , £ := mrlH , 

and set 

J C := C * = (curl it) o * = ( A T V) x (110$) = (A T V) xv, 

|£ := £0 * = (curl.ff) o * = (A T V) x (H o *) = (A T V) x B. 

Using the definition of the matrix ^4 in (7), the relations (36) can be easily inverted to find 

x d 3 v Vip x 5 3 J5 
curl u = C H j , curl B = £ H j . (37) 

Applying the curl operator to the original equations (2) satisfied by (it, H), we easily find that the Eulerian 
curls (C,£) solve the system 

f dtO* + (u* • V)C± - (i/ ± • V)|± - (C± • V)it± + (f± • V)F± = , 

ja*^ + (u ± • V)^ - (H+ • V)^ + [curl; it ± ■ V]iJ ± - [curl; H± ■ V]u ± = , 

in (J {t} x r2 ± (i). Making use of (36) and recalling the definitions in (9), it follows that (C, £) solve 
te[o,x] 

fc^ + (5* • V)C± - (S ± • V)£± - (AC± ■ V)«± + {At* ■ V)B± - , 

{5 t e ± + (^ ± -V)e ± -(i? ± -V)C ± + [ J 4 T Vx;Ai; ± -V]B ± -[ J 4 T Vx ; AB ± -V]i; ± =0, 1 J 
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in [0,T] x fi ± . Thus, in order to estimate the curl of d ± ,-B ± , we are reduced, after (37), to proving 
suitable bounds for the H 2 — norm of the solution (CO to (38). Let us observe that with our regularity 
assumptions on the original solution, there holds (CO ^(H 2 ) n C €(H Z ) so all integration by parts 
below arc legitimate. 

Let us introduce an associated energy functional Jfc defined by 

^(*) : =JEE/ {im^ + I^WI 2 }^. (39) 

1 ± \(>\<2 Jn± 

Differentiating with respect to t and making use of (7), (9), (38) gives 

■*"(*) = E E / {^ 5 *C ± ■ a^ + d^dtt ■ a ^} dx = jei(t) + j&{t) + jfr 3 (t) , (40) 

± | f 3|<2"' n± 

where 



J^(t) := - E E / { ( 5± ■ V ) 9 ^ ± - ( S± ■ V ) 9/3 ^ } ' ^ 

± |/3|<2^ 0± 

+ { (« ± • VJ^f* - (S ± • V^C* } • d ^ dx , 

^):=-EE / {[^;5 ± -V]C ± -[^;B ± -V]f t }.^c ± 

± |/3|<2^ 0± 

+ { [a' 3 ; 5* • V]^ - S ± ■ V]C ± } ■ 9' 3 e ± dx . 
-EE / ^ ((^ ■ V)5± - (^ ■ V)^) • 



± l/3|<2- 



([A T Vx;Aw ± • V] B ± - [A T Vx;AB ± • V] v ± ) • a^ dx . 



Let us estimate separately each of the above J^i, for i = 1, 2, 3. We start with To estimate this term, 
we use Leibniz' rule and integrate by parts. The boundary conditions (10) give 



J6(t) = ~E)£ / \t ± -v(\d\ ± \ 2 + \d^ ± \ 2 )-2B ± -v(d^ ± -dP( ± )}dx 

2 ± |0|<2 , ' n± 

= EE / (| div£i± (l^r + l^l 2 ) -div^a^-a^j dx. 



± |/3|<2" 

Applying Cauchy-Schwarz inequality, we obtain 

Vie [0,1"], |Jft(t)| <-F(<?(t)). (41) 

Let us now deal with the term J^. We focus on the first integral involved in the definition of J^2, 
namely 



E / [^;u ± -V]C ± -^C ± da;. 



In the sequel d l and d 2 stand for any derivative of order one and order two respectively. The commutator 
is zero if /3 = 0. If \/3\ = 1, the integral is of the form 

d 1 ^ d 1 v ± d 1 C ± dx. 

n± 

Using an L°° bound for d 1 ^ and Cauchy-Schwarz for the two remaining terms, we have 

lV± al„~,± QlA± 



a 1 ^ a 1 ?}* a 1 ^ dx 



<W)). 
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It remains to examine the terms in the commutator with \/3\ = 2. We can easily check that such a 
commutator can be rewritten as a sum of the form (we omit the harmless numerical constants) 

5 1 5 ± 5 2 C ± 9 2 C ± + d 2 v± <9 1 C ± <9 2 C ± dx . 

The first term is estimated as in the case |/3| = 1 by using an L°° bound for 9 1, ± . The second of these 
two terms requires more attention. We combine Holder's inequality and the Sobolev imbedding Theorem 
(recall that in three space dimensions H 1 is imbedded in L 6 ): 

Q2~± gl C ± d 2 ( ± dx < | a 2~±| 3 ^1^1 || 9 2 C ±|| ± < c || 5 ±|| || C ±||2 < F(<r(t)) _ 

n± 

In a completely similar way, we can handle the other commutators in J^(t) to finally get the estimate 

Vie[0,T'], |jfc(t)| < F(<f (t)) . (42) 
We now turn to the last term that we write in the form J^(t) = J^si(t) + Jfazif) with 

4(i):=-EE / ^ ' V ) B± - (^C ± • V)^) • dx , 

± |0|<2"' n± 

^32(t) ■= H f d f3 {[A T Vx;Av ± ■V]B ± -[A T Vx;AB ± -Vjv^ -d^^dx. 

± |/3|<2^ 0± 

The first integral in J%}i(t) are estimated by Cauchy-Schwarz inequality and by using the fact that 
H 2 (ft ± ) is an algebra: 



<||C ± lla I ±IH| 2 ||f t || 2 , ± ||B ± || 3 , ± <^(*)). 

The second integral in J^i (t) is estimated in the same way and we get 

Vte[0,T'], |jr 31 (i)| <F{g{t)). (43) 

As for J%} 2 (t), it is rather easy to sec that the quantity [i T Vx; Av ± ■ S/]B ± - [A T Vx; AB ± ■ V]^ can 
be expanded as a sum of terms of the form 

A d 1 A v ± d 1 B ± + A d l A B ± dh^ + AA d 1 v ± d 1 B ± , 

where we have disregarded the indices for the sake of simplicity. Hence the H 2 norm of this quantity can 
be estimated by a quantity of the form F(S'(t)). Using Cauchy-Schwarz inequality in J%32(t), we end up 
with 

Vie [0,2"], |JT 32 (i)| <iW)). 

Combining the latter estimate with (41), (42) and (43), we have obtained 

Vte[0,T'], \Jtr'(t)\<F(£(t)). 

We can now integrate this inequality from to t and use (37). The "error" terms Vip x d^vf/J, 
are estimated as in the paragraph on the divergence estimate, see (35), so eventually we 

get 

Vte[0,T'], || ( nirli; ± (f),(mlB ± (t)||l i± <Mo 

+ tF( max S («)) + C \\f(t)\\ 2 H2 , (T2) ^ (t) , B ± (t)\\l ± . (44) 

0<s<t K 1 
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3.3. Final estimate for the velocity and magnetic field. With the above divergence and curl 
estimates, we are ready to obtain the main a priori estimate for the velocity and magnetic field in each 
domain {l^ . The only point is to observe, through elementary algebraic manipulations, that the H 3 norm 
of a vector field is controlled by the L 2 norms of tangential derivatives of order < 3 and by the H 2 norms 
of its divergence and of its curl. We thus add the estimates (34), (35) and (44) to obtain 

Vte[0,T'], \\v ± (t),B ± (t)\\l± <M +tF(max/( S )) + C ||/^^ 

0<s<t v ' 

where, of course, the numerical constants Mo, Co arc independent of the solution. Consequently, up to 
choosing Eq small enough so that Co £o < 1/2 an d adapting the time interval [0, 2"] so that (17a) is valid 
with the new definition of £o> we obtain 

Vie [0,2"], \\v ± (t),B ± (t)\\l ± < M +tf(max S{s)) . (45) 



4. Estimate of the front 
From the linear system of the boundary conditions on T 

(B+d 1 f + B+d 2 f = B£, 
\b 1 -9 1 / + j B 2 -9 2 / = B 3 -, 

we have already seen that the determinant B^ B 2 — B 2 B^ does not vanish on [0, T'] x T. More precisely, 
we have 

\ Bl B 2 -B 2 B x (t,z,0)\ - i + \ V >f(t, X 'W -4(1 + Ce 2 )' 

where we have used (18a), (17a) and the imbedding 7? 15 (T 2 ) L°°(T 2 ). We also note that thanks to 
(17b), the L°° norm of B^ is uniformly controlled on [0, T']. Therefore, using the latter uniform bound 
for the determinant and inverting the linear system (46), we have 

Vie [0,2'], ||V7Wlk^(T=) < Co ||B ± (*)||3,± , (47) 

with Co depending only on So and R. 

From the other boundary conditions on T: 

dtf = vf - vf dif - vf d 2 f , 
(47) and the fact that 72 2 - 5 (T 2 ) is an algebra, we infer the second main estimate for /: 

Vie [0,2'], ||a t /(*)|| fla .» (r i) < Co (||« ± (*)||3,± + ||« ± (*),B ± (t)||i,±) ■ (48) 
In particular, we can integrate from to i and get 

Vie [0,2"], \\f(t)\\ H 2,, iJ 2 ) <\\fo\\ H 2, HJ2) +tF(m^x^(s)). (49) 

0<s<t 

We simplify (47), (48) and (49) by using (45) (we feel free to use t 2 < t which always holds by assuming, 
without loss of generality T" < 1): 



Vie [0,2"], \\d t f(t)\\ 2 H2 . 5(r2) <M +tF(max £{s)) , 



0<s<t 

\\f(t)\\ 2 H3 . 5{T2) <M +tF(m & x #(s)), (50) 
11/(4)11^.5^2) < \\f \\ H 2. 5(j2) +tF(max t £(s)). 

The last estimate in (50) says that /(t) remains small in H 2 5 provided that we start from small initial 
data and the first and second estimates in (50) give a control of dtf(t) in H 2 5 and / in H 35 . We observe 
that fit) is expected to remain small in H 2 5 but has no reason to be small in H 3 5 (in particular because 
no smallness condition has been made on the norm of fo in H 35 ). 
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5. The elliptic problem for the total pressure 

We first deduce from (8) the elliptic system of equations solved by the total pressure. Applying A T S7- 
to the equation for in (8) gives 

-A T V ■ {A T V Q ± ) = A T V • {c> t w ± + (« ± • V)« ± - (B ± • V)B ± } . 

Using the divergence relations A T \7 ■ = A T \7 ■ B ± = 0, we then deduce the equations 

— A T V ■ (A T V Q^) = , (51) 

where we have set 

JT± := -d t A kl d k vf + A kl d k v± ■ Vvf - 5 ± • VA kl d k vf - A kl 9 fc B ± • VBf + B ± • VA kl d k Bf . (52) 

Recalling that a = J A we get from (51) the equivalent equations 

-a T V ■{A T VQ ± ) = J&±. (53) 

Now we look for the boundary conditions satisfied by Q ± . Since = Bf = and = vf = Bf = on 
[0, T] x T±, from the third equation for in (8) evaluated on we obtain the homogeneous Neumann 
condition 

d 3 Q ± = on[0,T]xr±. (54) 
On r we take the scalar product of the equation for v in (8) with the vector N. We get 

- (A T VQ ± ) • N = {dtv* + (S* ■ V)v ± - (B* ■ V)S ± } ■ N . (55) 

Let us compute the jump of each quantity in (55) across T. Since [Q] = gives [d\Q] = [c^Q] = on 
[0, T] x T, we obtain (recall that J = 1 on T) 

[(A T VQ) ■ N] = [A tj Nj d t Q] = (1 + |V'/| 2 ) [d 3 Q] . (56) 

Using the boundary conditions dtf = v ± ■ N, B ± ■ N = 0, on [0, T] x T, we also deduce 

{d t v + (v ■ V)v - (B ■ W)B} -n] = [2v'- Vd t f + (v' • V')V/ • v' - (B' ■ V')V/ • B'] . (57) 

Thus from (55), (56) and (57), we find the boundary condition 

[A tj Nj diQ] =<S on [0, T] x T , (58) 

where we have set 

Sf := - [2v' ■ V'd t f + (v' ■ V')V'/ ■ v' - {B 1 • V')V/ ■ B'\ . (59) 
Collecting the equations (51), (54), (58) gives the elliptic problem 

'-A T V- (A T VQ ± ) = jr± t on[0,T]xO ± , 

[Q] = 0, on[0,T]xr, 

[Atj Nj d e Q] = &, on [0, T] x T , (60) 

d 3 Q ± = on[0,T]xr±, 

(xi,X2) i-> Q ± (t, Xi,X2,Xs) is 1 — periodic, 

with ^ ± and 'S defined in (52), (59), respectively. 

Remark 6. When one tries to solve the elliptic system for the pressure, it may be easier to work with the 
formulation (53) instead of (51) because of the necessary compatibility condition on the data t ^" ± , ^ '. More 
precisely, trying to solve problem (8) by a fixed point argument, one possible step could be the resolution of 
system (60). (We have in mind the approach used in [11], for the resolution of the incompressible MHD 
equations in a fixed domain under slip boundary conditions.) Thus the compatibilty condition needs to be 
satisfied by the data. 
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In order to formulate the compatibility condition we compute by an integration by parts 

-V / a T V • (A T VQ ± )d.T = - / a 3i A ki d k Q + Ax' + / a 3i A ki [d k Q] Ax' + / a 3l A kl d k Q~ Ax' 
± Jn± Jr + Jr Jr_ 

+ / d k a ki A hl d h Q ± Ax, 

where the last integral vanishes because of the so-called Piola's identity d k a k i = 0. The boundary condi- 
tions for Q yield 

- V / a T V • (A T VQ ± ) Ax = [ a 3l A kl [d k Q] Ax' = [ A ki N t [d k Q] Ax' . 
± Jn± Jr Jr 

This shows that the data & , ^ of problem (60) need to satisfy the condition 

V / J& ± Ax= f Sfdz'. 
^7n± Jr 

This condition is satisfied with our definitions since 

V / J.^ ± Ax = Y ( a T V • {d^ + (u ± ■ V)u ± - (5 ± • V)B ± } Ax 
± Jn± ± Jn± 

= - J N ■ {d t v + (v ■ V)v - (B ■ V)£?} Ax' = J&Ax', 

from (57), (59), and by computations as above. Thus the compatibility condition is satisfied. 

Our approach here is different because we have already assumed that the solution exists and we only 
wish to prove an a priori estimate on a time interval that is independent of the solution. Consequently, 
we shall deal with the slightly more symmetric formulation (51) to derive energy estimates. 

In the rest of this section we study the elliptic problem (60) for generic data J^^Sf. Only at the 
end of the section we will go back to the specific definition of J^,^ given in (52), (59). As (60) is 
time-independent, in the sense that time appears only as a parameter, for simplicity of notation from 
now on in this section the explicit dependence on t will be neglected. 

5.1. The functional framework. Thanks to the continuity of the total pressure across T, we can define 
the pressure Q € i/ 1 (0) by Q := on il ± . The function Q belongs to the Hilbert space 



R e H 1 ^) , J RAx = oj 



The space "V equipped with the norm ||V.R||l2(q) is indeed a Hilbert space, because of the Poincare 
inequality, and the norm || Vi2||£2(n) is equivalent to the standard H 1 norm. In what follows, the function 
Q will be estimated in the space and we shall repeatedly use the fact that the L 2 norm of VQ is 
equivalent to ||Q ||i,±- 

5.2. The general procedure for the pressure estimate. Step 1 We start from (60), multiply each 
equation in by Q^, integrate over il 1 * 1 and use integration by parts. This yields 

J2 I d k (A kj Q ± )A ej d e Q ± Ax= [ A 3j Q+ A ej d e Q+ Ax' - [ A 3j Q~ A ej d e Q~ Ax' 
± Jn± Jr + Jr_ 

-J A 3] Q+ A ei d i Q+Ax' + J A 3j Q~ A t , d t Q~ Ax' 
&± da- . 
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We recall that from the boundary conditions, ip and d^Q vanish on T± so the integrals on T± vanish. 
So we get 

Y [ A k j duQ* Mi deQ ± dx = - V / {d k A kj ) Q ± A (j d e Q ± dx 

- / Qlr^dx' + V / Q ± ,? ± dx, 
Jr ± Jn± 

where Q|r denotes the common trace of Q on T. The integral on the left hand side gives the coercive 
term in VQ ± (see the definition (7) and recall the condition ||VV , ||LooQo,T']xf2) < 1/2)- Then we apply 
the Cauchy-Schwarz and Poincare inequalities to derive 



c\\Q 



+ |2 



: \\» 



«-±l|2 



^iiff-o. 5(T =)+E / \ Q * A ^\ 10*1 1^1 dx, 

for a suitable numerical constant c > 0. Then we use the Holder and Sobolev inequalities to derive 
V / |9 fe ^||Q ± ||^Q ± |dx<C||VQ ± || ± |VA| 4 |Q ±l 



4,± 



<CP|| 2 ||Q ± ||? ± <C||/(t)|| ff2 .5 (T2) ||Q ± ||? i± , 



Up to choosing eg small enough, we have thus derived the first estimate 



Vt e [0,T'] , ||Q±||2 ± < Co ( H^lli + HSfll^-o..^) ) ■ 



(61) 



Step 2 We are now going to estimate Q in H (f2 ). Let us first apply a tangential derivative 9 to 
(60), with d = d\ or d = d 2 - Defining cf 1 := 9Q ± , we obtain the elliptic system 

' -A T V ■ (A T VQ^ = W ± , on [0, T] x fi± , 



[Q] = , 
fcQ* = 



on [0, T] x T , 
on [0,T] x T, 
on [0,T] x r± , 
(xi,x 2 )^Q {t,xi,x 2 ,x 3 ) is 1 - periodic, 

where the new source terms t ^" ± , Sf are defined by 

^ - + dA kj d k {A tj $Q±) + A kj d k ((EA ej ) fifcQ*) , 
¥:=dif- d(A ej Nj) [diQ] =m- d(\Vf\ 2 ) [d 3 Q] . 
We apply the same procedure of integration by parts as above, obtaining first 

Y [ A k3 9 fe Q ± Mi dx = - V / (d k A kj ) Q ± Ay fl^ dx 
± ./n± ± 7n± 

- / Q|r^dx' + V / Q^^dx. 



(62) 



(63) 
(64) 



where Q|r denotes the common trace of Q on T. The integrals on the left hand side give the coercive 
terms and, as above, we can absorb the first integrals occuring in the right hand side by choosing £0 small 
enough. We thus have 



c||g ± ||i± < - / Qlr^dx' + V / Q ± ^ ± dx. 
Jr _i_ Jn± 



We now estimate the integrals on 17 ± , recalling the definition (63) for . Let us first observe that 
the term with d^ ± can be integrated by parts and we can then apply Cauchy-Schwarz and Young 
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inequalities. The other terms are estimated as follows: 

Y [ IQ* | \dA kj I \Atj | \d k d t Q ± | dx^CWQ* || 2 ,± |VA| 4 |A|oo IQ^i 

<C||A||l||Q ± ||l i± <C||/(t)||| P .» (Ta) ||Q ± ||l i± , 

V / \Q ± \ \dA kj \ \d k A ej \ \d t Q±\6x <C^t\^ ± \VA\l |Vg ± | 4 ,± 

< C mill ||Q ± ||l i± < c ||/(t)||^. 5(T2) ll<5 ± lli,± , 

and applying similar sequences of inequalities, the reader can get quickly convinced that all other terms 
in the product Q & are estimated by the same quantity. We thus have 



c||Q ± ||? !± <||^ ± || 2 ± 



Qlr^dx' 



C||/(t)||| ra . B(Ta) ||Q 



±||2 



I2,± • 



Let us now turn to the boundary term. Of course, we have 



Q\ r d^dx' 



< 



||jJ0.B(T2) ||Q|r||jJO.B(r) < C \ft\nQ-*<$*) ||<3 || 1)± . 



The remaining term occuring in is easily estimated as follows: 



Q\r[dzQ]d{\V'f\ 2 )&x' 



< \Q\rl \[d 3 Qh \d(\Vf\ 2 )\ 



<C||Q| r ||Ho. 5(r) ||[^Q]|| ff o. 5(r) |||V7rllifi-* m ^IIQ^Ilill/WII^"^)' 
where we have used H°- 5 (T) ^ L 4 (r) (which holds in two space dimensions), and the fact that i? 15 (r) 
is an algebra. Applying Young's inequality again, we thus obtain 

\\Tf\\l ± < Co (W^Wl + \\n 2 Ho. 5(T2) + \\f(t)\\ 2 H2 . HJ2} ||Q ± ||5 i± ) . (65) 

Step 3 The remaining second order derivative d 2 Q ± is estimated directly from the equation (60) by 
using the explicit expression of the coefficients A k j- More precisely, (60) reads 

Aji A kl fydkQ* = - Ajt djA u , 

that is, 

1 + |VVf 2 , 2 , 2 . dijjdid^ n d 2 jj d 2 d 3 Q ± _ ± ± 

(TTW sQ l0 2 ° " 2 i + 9 3 ^ " i + • l '" vu '' y,g • (66) 

We thus obtain 

cii^nl < c (ll^llU + H^Hi + HAS^a^Hi 

< c (iioliU + n^iil + ||/(i)ll^. 5(T =) ||Q ± ||5 )± 

Combining with (61) and (65) and choosing the numerical constant Eq sufficiently small, we obtain 

Vie[0,T'], \\Q^ ± < Co (11^111 + HSfHlo..^)) ■ (67) 

Step 4 We now apply the estimate (67) to the solution Q ± to the problem (62), which has the same 
form as (60) but with different source terms (defined in (63) and (64)). We thus have 



— ±„9 , ~ /,— ±„_ 

ff 0.5( T 2) 



Vte[0,T'], \\CfM ± <c{\\& \\l + \\n 



The L 2 -estimate of follows by applying similar arguments as above; for instance, we have 
||aWA9 1 Q+|| + < Wd'Ad'AU ||Q + || w i.- (n+) < C |VA|1 ||Q+|| 3>+ < C \\f(t)\\U*<r) ll^lk 
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All the other terms in admit the same upper bound, that is 

HJ^Hi < c + c ||/(t)||^. 5(T2) ||Q ± ||i 1± ) . 

As far as the boundary source term is concerned, we apply Lemma 5 and obtain 

P(|V7| 2 ) [d 3 Q]\\ H o, ( r) < C p(|V7| 2 )lk-( T s) \\[dsQ}\\ H i. H r) < C \\f(t)f H ™m l|Q ± lli 1 ± ■ 
We have thus derived the upper bound 

vt g [o,t'] , iiq*!!^ < c (n^im + wn 2 H^ m + \\f(t)\\'U H T*) wqHu) ■ 

The remaining third order derivative d^Q^ can be estimated by applying ^3 to the equation (66). The 
commutators are estimated exactly as above, and we now feel free to skip a few details. Eventually, up 
to choosing a sufficiently small numerical constant £0 > 0, and provided that 2" is such that (17a) holds, 
we derive the estimate 

Vi G [0,2*] , HQ^l'i < Co (ll^HU + ||Sf|||ri..(T»)) ■ (68) 

5.3. The final pressure estimate. It only remains to use the definition of the source terms , ^ in 
(68). Using first the fact that i2 1-5 (T 2 ) is an algebra and recalling the definition (59) of Sf , we have 

\\&(t)\\ m , m < C (||« ± (t)||3,±||9 t /(t)llff'-«(P) + \\v ± (t),B ± (t)\\l±\\f(t)\\H^)) , 
and using (45), (50), we get 

\W(t)\\ 2 Hl _ HT2) < M + tF(max <?(*)) . 
The source terms can be estimated by applying the classical estimate 

IK «2||#i < C (|K||l=° \ W2WH 1 + IKIU=° II ff i ) • 

Analyzing each separate term in the definition (52) of by applying the latter product estimate and 
by using (17), (45) or (50), we get 

11^ Wlli± < M + tF(max § (s)) . 

0<s<t 

Adding the previous two inequalities, we obtain our final estimate for the pressure: 

VtG[0,2"], ||Q ± ||| ± < M +tF( max #(s)) . (69) 

0<s<t 

6. Proof of Theorem 4 

If we summarize the analysis of the previous sections, we have shown that there exist some numerical 
constants £0 > and Mo > 0, there exists a nonnegative nondecreasing function F on R + , all three 
depending only on Sq and R such that, on any time interval [0,2"] for which the inequalities (17) are 
valid, there holds 

Vi G [0, T'\ , £{t) <M + t F( max £{s)) . (70) 

The function F and the constants £0, M are independent of the particular solution that we are consid- 
ering. Moreover, i2 2 (Jl ± ) is an algebra so applying direct estimates on (8) we find 

yte[o,r], iift« ± (t),a t B ± (t)i| 3l ±<i ? (^(t)) > 

so integrating with respect to t we have 

Vte[0,T'], \\v ± (t)-v t ,B ± (t)-B t \\ 2 ,±<tF{max £(s)). (71) 

0<s<t 

From now on, the nonnegative nondecreasing function F is fixed, as well as the constants £0, Mq. 
To complete the proof of Theorem 4, we define E\ := Eq/2, and we choose a time To > such that 
2 T F(2 Mo) < M and 2T F(2M ) < £1. We emphasize that the definition of T only depends on S 
and R. Then we define 2" as the maximal time on which (17) holds (2" is positive because (17) holds at 
the initial time with a strict inequality). We will see that T < 2" if T < T, and 2" = T < T if T < T . 
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There arc now two possibilities. Let us first assume T > T , and let us define / as the set of all times 
t G [0, T ] such that 

max S'(s) <2M , max ||w ± (s) - v k 7 B ± (s) - B^\\ 2 .± < e , max \\f(s)\\ H 2.5„2 ) < e ■ 

0<s<t 0<s<£ 0<s<t 

Then / is non-empty since it contains (use (70) for t = 0), and I is closed since all functions involved 
in the definition of / are continuous. Let us show that I is open. Let t£ I. Using (70), we have 

<g(t) < M + tF( max g(s)) < M + T F(2 M ) < 2 M . 

~ 0<s<* 

In the same way, (50), (71) and the definition of e\ give 

||» ± (t) -v k ,B ± (t)-B±\\ 2 ,±<e , \\f(t)\\ H ^(^) <e - 

Consequently, there exists a neighborhood of t in [0, To] that is included in I. In other words, / is open. 
Hence I = [0, To] and the result of Theorem 4 is proved. The proof in the case T < To is similar. 



7. Proof of Lemma 1 

Given x G Cg°(K), x = 1 on [-1, 1], we define 

fW(x',x 3 ):= X (x 3 \D\)f(x'), i>(x',x 3 ):=(l-xl)fW(x',x 3 ), (72) 

where ^(x3|D|) is the pseudo-differential operator with \D\ being the Fourier multiplier in the variables 
x' . From the definition it readily follows that ip(x',0) = f(x'), ^>(x',±l) = for all x' G T 2 . Moreover, 

d 3 ^(x',x 3 ) = -2x 3 f^\x',x 3 ) + (1 - x 2 3 ) X '(x 3 \D\) \D\ f(x') , (73) 

which vanishes if x 3 = 0. Given any function g defined on T 2 , let us denote by Ck{g) the fc-th Fourier 
coefficient 

c k (g)= f e~ 2ivk - x ' g(x')dx', fceZ 2 . 
Since c fc (/ (1) (-, x 3 )) = x{x 3 |fc|)c fc (/), we compute 

M;x 3 )\\ 2 Hm(T2) = (1 - x 2 ) 2 \\f W (-,x 3 )\\ 2 Hm(T2) <C(1- x\f Y, (1 + l fc H 

fcez 2 

< c (i - x\f J2 c 1 + i fc i 2 r * 2 (*3 i*d m/)i 2 ■ 

fcez 2 

It follows that 

Wij 3(ff »(x-))<C f {l-xlf Y( 1 + \ k ft m X 2 (x 3 \k\)\c k (f)\ 2 dx 3 
<C ^(l + |fc| 2 r|c fc (/)| 2 f X 2 (x 3 \k\)dx 3 

<C|c (/)| 2 +C £(l + |fc| 2 r| Cfc (/)| 2 -L r x 2( s )ds. 

\k\>l ' ' ^~l fe l 

Denoting by I £ C°°(R) the primitive function of x 2 vanishing at — oo, i.e. X'(s) = X 2 {s), we notice 
that X is bounded over all R. Then 

HWh ( ^(T 2 )) < C |co(/)| 2 + C ]T (1 + |fe| 2 r- 1/2 M/)| 2 sup |X(s)| < C \\f\\ 2 Hm -^ m ■ (74) 

lfcl>l sm 



C k (fV(;X 3 )) 
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In a similar way, from (73), we obtain 

IIMIi^-ifr*)) < G (llx(*3 I^D /lli;,^-^)) + Hx'fe \B\) \D\ /||ij 3(ifm _ 1(T2)) ) 



i 



+ C £(l + |fc| 2 r-M*| a M/)| a T Ix'^slfeDl^a 

<cn/n| w _, /aCP) +cx;(i+i*i a r- 1 i*iicfc(/)i 2 /' fel ix( S )i 2 d S . 



Denoting by y G C°°(IR) a primitive function of (x') 2 , wc also notice that Y is bounded over all K, so as 
in (74), we get 

||9 3 V||| 2 ( ^- 1(T =)) < C ||/||| m -, /acP ) + C ^ (1 + Ifcl 2 )™" 1 / 2 |c fc (/)| 2 sup |Y( S )| < C \\.f\\ 2 Hm - l/2{J2) . 

\k\>i sm 

Iterating the same argument yields 

W^Wll^H^-i^)) ^ C|i/llff m -i/2 (T 2) , j = 0,...,m. 

Adding over j = 0, . . . , to finally gives ?/> € H m (Q) and the continuity of the map / i— > -0. 

The proof of Lemma 2 follows from Lemma 1, with t as a parameter. Notice also that the map 
/ — > see (72), is linear and that the time regularity is conserved because, with obvious notation, 
(djfY 1 ^ = d^(f^). The conclusions of Lemma 2 follow directly. 
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